>1 

AD-A023  231 

STEADY  STATE  ANALYSIS  OF  HYDRAZINE  CATALYTIC 
THRUSTERS  FOR  DIFFERENT  TYPES  OF  CATALYSTS 

Antonio  Crespo 

Instituto  Nacional  de  Tecnica  Aeroespacial 


i 


Prepared  for: 

Air  Force  Office  of  Scientific  Research 
15  January  1976 


DISTRIBUTED  BY: 

Knn 

National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 


I 


1 


SECURITY  CLASSIFICATION  C*F  THISJJACE  fHTi.o  O at*  F.ntrted) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER  ~ - ... 

AFOSR  - TR  - 7 0 - 013  1) 


D A ftp  READ  INSTRUCTIONS 

, t-  - ‘ - BEFORE  COMPLETING  FORM 

2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMUCH 


4.  TITLE  (arid  Subtitle)  * 

STEADY  STATE  ANALYSIS  OF  HYDRAZINE  CATALYTIC 
THRUSTERS  FOR  DIFFERENT  TYPES  OF  CATALYSTS 

7.  author^ 

ANTONIO  CRESPO 


5.  TYRE  OF  REPORT  A PERIOD  COVERED 

INTERIM 

1 Jan  1975-30  June  1975 

«.  PERFORMING  ORG.  REPORT  NUMBER 

Report  #4 

4.  CONTRACT  OR. GRANT  NUMBERfx) 

AFOSR  72-2253 


>.  PERFORMING  ORGANIZATION  NAME  AND  ADORESS 

INSTITUTO  NACIONAL  DE  TECNICA  AEROESPACIAL  (INTA) 
PASEO  PINTOR  ROSALES  34 
MADRID  8 SPAIN 

11.  CONTROLLING  OFFICE  NAME  AND  ADORESS 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NA 
BUILDING  410 

BOLLING  AIR  FORCE  BASE.  D C 20332 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 

681308 

9711-01 

61102F . 

12.  REPORT  DATE 

Jan  1976 

13.  NUMBER  OF  PAGES 

yt 


147'  MONITORING  AGENCY  NAME  4 AOORESSf It  dlllerent  from  Controlllni  Olllc e)  IS.  SECURITY  CLASS.  ( ol  thl,  report) 

UNCLASSIFIED 


ISn.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


|16.  DISTRIBUTION  STATEMENT  (ot  thlo  Report) 


Approved  for  public  release;  distribution  unlimited. 


J *7.  DISTRIBUTION  STATEMENT  (ot  tht i obetroct  entered  In  Block  20,  It  dlllerent  from  Report )' 


14.  SUPPLEMENTARY  NOTES 


19*  KEY  WOROS  (Continue  on  rtvaraa  aide  It  nacaaaary  and  Identity  by  block  number) 

HYDPAZINE  THRUSTERS 
CATALYSTS 

AMMONIA  DECOMPOSITION 

ASYMPTOTIC  ANALYSIS  0 


PRICES  SUBJECT  TO  CHANGE 


20.  ABSTRACT  (Continue  on  rarer**  aide  It  nacaaaary  and  Identity  by  block  number) 

A theoretical  analysis  has  been  made  to  determine  the  characteristics  of  the 
steady  state  behaviour  of  a hydrazine  decomposition  engine.  The  model  makes  use 
of  the  simplification  that  consists  in  considering  that  the  ammonia  decompositio 
is  only  important  at  high  temperatures. ' The  results  are  valid  for  different 
types  of  catalysts  with  different  stoichiometries  and  reaction  rates  for 
hydrazine  and  ammonia  decomposition.  Conditions  at  the  chamber  exit  are  obtaine 
analytically  as  a function  of  chamber  length  and  other  characteristics,  in 
particular  expressions  are  given  for  fractional  ammonia  decomposition  at  the 


FORM 
I JAN  73 


EDITION  OF  1 NOV  45  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wien  Data  Entered) 


HTHeEtSSIOTr 


tCCUBlTV  CLASSIFICATION  OF  THIS  PAtf  CfHTitn  £>«<•  Knit  ltd) 


f ^ . A** 

l exit  and  the  specific  impulse  obtained  if  the  outgoing  mixture  would  be  expanded 
to  the  vacuum.  Fluid  and  particle  temperature  distributions  along  the  chamber 
are  also  calculated  keeping,  however,  fixed  a nondimens ional  parameter 
experssing  the  ratio  of  the  mass  and  heat  transfer  from  the  interstitial  fluid 
to  the  particle. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEWicn  Oil*  Enfttd) 


£*“*««y*' i**»1 


1 13084  j> 

AFOSR  ” TfU  7 6-0186 


Giant  Humbe.fi.  AFOSR  72-2253 


STEADY  STATE  ANALYSIS  OF  HYDRAZINE  CATALYTIC  THRUSTERS 
FOR  DIFFERENT  TYPES  OF  CATALYSTS  • 


A.  Cl e&po 

Institute  Nacional  de  Td.uu.ca  Ae.fi.oc.4pact.al 
"E4te.ban  Te.fi.fi.ada4",  Madfiid,  Spain. 


15  Januafi.it  1976 


InteAim  Scientific  Repo At  *4 

1 Januaiy  7975  to  If  June  19  75 


Piepaied  foA 


D D Cv 

TpEEDinEIj 
Uj| \?R  2Q  1976 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  - M 

A Alington  Viiginia 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  (A?SC) 
NOTICE  OF  T5WTCKII2AL  TO  DDC 

, This  tech.  leal  i * lux*  been  reviewed  and  is 

and  approved  ;'*•»*  . -1  c r:  lease  IA./  A FR  190-12  (7b). 

Ditiii’iL'Utio.i  ia  ualiraited. 

A«  D. 

Xoohnical  Information  'Officer 

EUROPEAN  OFFICE  OF  AEROSPACE  RESEARCH  A HV  DEVELOPMENT 
London,  England. 

, " KMnoutcO  »y 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U.  S.  OEMRTMENT  OF  COMMERCE 
SPRIMFIElp,  VA.  22111 

Approved  for  public  release ; 
distribution  unlimited. 


V 


t 


ABSTRACT 


- / A theoretical  analysis  has  been  made  to  determine  the 
characteristics  of  the  steady  state  behaviour  of  a hydrazine 
decomposition  engine.  The  model  makes  use  of  the  simplification 
that  consists  in  considering  that  the  ammonia  decomposition  is 
only  important  at  high  temperatures.  The  results  are  valid  for 
different  types  of  catalysts  with  different  stoichiometries  and 
reaction  rates  for  hydrazine  and  ammonia  decomposition.  Conditions 
at  the  chamber  exit  are  obtained  analytically  as  a function  of 
chamber  length  and  other  chatacteristics , in  particular  expressions 
are  given  for  fractional  ammonia  decomposition  at  the  exit  and  the 
specific  impulse  obtained  if  the  outgoing  mixture  would  be  expanded 
to  the  vacuum.  Fluid  and  particle  tempo. ature  distributions  along 
the  chamber  are  also  calcula+ed  keeping,  however,  fixed  a non- 
dimensional  parameter  expressing  the  ratio  of  the  mass  and  heat 
transfer  from  the  interstitial  fluid  to  the  particle. 
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INTRODUCTION  j 

The  purpose  of  this  work  is  to  carry  out  a theoretical  I 

analysis  that  will  provide  a set  of  relations  which  can  be 
used  as  a basis  in  a first  step  of  the  design  0f  hydrazine 
engines  with  catalytic  decomposition.  In  this  work  we  shall 
only  consider  the  steady  state  behaviour  of  the  engines. 

Previous  work  on  the  theoretical  aspects  of  this  problem  has 

1 2 

been  made  by  Kesten  and  our  research  group  of  INTA  . Works 

of  a more  practical  nature, including  experimental  results,  have 

3 4 

been  made  by  the  Rocket  Research  Corporation  and  Adler  et  al  . 

. 

| 

In  all  these  works  and  in  the  present  one,  it  is  assumed 
that  Hydrazine  decomposes  yielding  Ammonia,  Hydrogen  and  Nitrogen, 
and  in  turn  Ammonia  decomposes  also  giving  Hydrogen  and  Nitrogen. 
Whereas  in  the  previous  works  it  was  assumed  that  the  stoichiometry 
of  Hydrazine  decomposition  was  given,  in  the  present  work  we 
consider  it  as  a parameter  making  the  results  more  general. 

Kesten  obtained  extensive  numerical  results  and  provided 
a computer  program  for  further  calculations  (Smith,  Smith  and 
Kesten5).  Although  he  also  gave  some  correlations  of  has  nume- 
rical results  for  some  range  of  variation  of  the  parameters,  it 
is  always  difficult  to  extend  these  correlations  to  other  para- 
meters of  interest  and  to  obtain  the  physical  insight  of  the 
probl . n that  an  analytical  solution  provides.  For  example  his 
correlation  is  only  valid  for  the  Shell  405  Catalyst  and  some 
particular  conditions  that  were  assumed  regarding  stoichiometry,- 

velocity  of  reaction  diffusion  coefficients  etc..  The  previous  I 

0 1 

work  made  at  INTA  made  some  simplifications  on  the  model  of  j 

Kesten  tha-c  did  not  affect  essentially  the  physics  of  the  problem,  I 

obtaining  in  turn  a solution  with  all  the  characteristics  of  the  | 

reactor,  that  depended  only  on  three  non-dimensional  parameters  j 


/ , 
of  which  only  ;oye  appeared  to  be  of  importance.  This  model 

still  retained  some  parameters  fixed  such  as  the  kinetics  of 
-hy/lrazine  decomposition-  and-  the  reaction  rate  of  ammonia  decom- 
position. The  present  model  is  a further  simplification  of  the 
2 

previous  work  that  essentially  consist  in  considering  that 
ammonia  decomposition  is  only  important  at  temperatures  close 
to  the  maximum  one  that  is  attained  in  the  reactor.  This  simpli- 
fication, in  turn  allows  to  retain  as  variable  parameters  of  the 
-problem  all  the  characteristics  determining  the  kinetics  inside 
the  particle  and  the  transport  coefficients,  making  the  results 
appliable  to  other  types  of  catalysts. 


In  the  first  part  of  the  work  the  problem  is  stated  and 
a review  is  made  of  the  formulation  of  the  problem  as  was  pre- 
sented in  ref.  2,  extending  it  in  order  to  include  more  general 
characteristics  of  the  catalyst.  In  the  second  part  the  problem 
is  solved  considering  that  there  is  no  ammonia  decomposition, 
providing  a first  approximation  to  some  aspects  of  the  solution 
of  the  problem.  In  the  third  part  the  ammonia  decomposition 
effect,  is  retained. 


The  temperature  and  concentration  distributions  along  the 
chamber  are  calculated  as  functions  of  the  particle  surface 
temperature  distribution,  and.  in  order  to  calculate  it,  it  is 
necessary  to  solve  a single  ordinary  differential  equation  of 
first  order  in  which  appear  two  parameters;  one  related  to  the 
ammonia  decomposition  rate  and  the  other  being  the  ratio  of  the 
heat  and  mass  transfer  coefficients  from  the  interstitial  fluid 
to  the  particle.  A solution  of  this  equation, appropiate  for 
chamber  exi.t  conditions,  is  found  that  gives  temperature,  concen- 
tration of  species, and  specific  impulse.  This  solution  reproduces 
most  of  the  aspects  of  Kesten's  correlation.  However,  in  order 
to  find  other  parameters  such  as  maximum  temperature  and  its 
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location  in  the  chamber  it  is  necessary  to  solve  that  differential 
equation.  The  nature  of  the  solution  of  this  differential 
equation  depends  critically  on  whether  the  parameter  representing 
the  ratio  of  the  mass  and  heat  transfer  coefficients  +0  the 
particle  is  smaller  or  larger  than  one;  in  ref.  2 it  was  assumed 
to  be -equal  to  one,  whereas  in  ref.  1 it  was  taken  to  be  less 
than  one.  The  most  fundamental  difference  in  the  results  seems 
to  be  that  for  this  number  smaller  than  one  a maximum  is  reached 
in  the  particle  surface  temperature,  whereas  if  it  is  equal  or 
largerthan  one, the  particle  surface  temperature  always  decreases. 
Since  the  most  important  results  for  the  designer  are  those 
corresponding  to  the  chamber  exit  , and  these  are  obtained  anali- 
tically,we  do  not  dwell  much  in  the  present  work  on  the  signifi- 
cance of  this  parameter,  and  reserve  its  discussion  for  future 
research.  Numerical  solutions  are  presented  for  this  parameter 
equal  to  0.6  that  is  approximately  the  value  given  in  ref.  1. 


1.  FORMULATION  OF  THE  PROBLEM 

The  problem  to  be  studied  consists  on  the  steady  state 
analysis  of  ^the  catalytic  chamber  of  a hydrazine  microthruster 
for  space  propulsion. 


In  the  chamber  hydrazine  flows  between  the  interstices  of 
the  particles  and  diffuses  to  their  interior, where  it  decomposes 
generating  heat  and  products  that  are  transferred  to  the  inters- 
titial fluid.  In*  the  kinetic  model  assumed  here  hydrazine 
decomposes  according  to  the  reaction 


2N2H4  2(2-a)NU3  + ctN2  +C3a  - 2)»2 


and  in  turn  ammonia  decomposes  also 


2HH3  - N2  + 3H2 


The  hydrazine  decomposition  process  is  exothermic  with  a heat 
of  reaction  per  unit  mass  of  reactant  H, 


Hh  = 1930  + 300(l-a)BTU/lb 


and  the  ammonia  decomposition  is  endothermic  with 


H = 1400  BTU/lb 

a 


for  temperatures  around  1100°K  appropiate  for  this  study. 


In  refs.  1 and  2 it  was  assumed  that  a = 1,  and  in 
ref  3 a=2/3  for  the  Shell  405.  In  ref.  4 a was  taken  2/3  for 
other  type  of  catalyst. 


presented  schematically  in  figure  1.  Three  regions  of  the 
microthruster/  are  considered.  In.  a first  vapor  region 
:hydrazine:  enters  the  chamber  in  liquid  form  and  decomposes 
slowly  keeping  its  feeding  temperature  until  it  suddenly 
increases  to  the  boiling  temperature  at  the  operating  pressure. 
Then,  follows  a region  where  liquid  hydrazine  and  its  vapors 
and  those  of  the  decomposition  products  coexist.  In  this 
region  hydrazine  is  being  vaporized  and  the  temperature  of 
the  mixture  does  not  change.  When  all  the  hj drazine  is  vaporized 
in  the  third  vapor  region,  the  heat  released  in  the  hydrazine 
decomposition  goes  into  heating  the  mixture  increasing  its 
temperature.  This  temperature  increase  continues  until  the 
endothermic  ammonia  decomposition  process  becomes  more  important 
than  the  exothermic  hydrazine  .decomposition  process,  after  that, 
when  there  is  almost  no  hydrazine  left,  the  temperature  of  the 


flow  decreases. 


In  the  previous  work  carried  out  at  INTA  the  following 
equations,  in  the  vapor  region,  expressing  that  the  variation 
rate  of  the  different  species  along  the  chamber  is  equal  to 
the  rate  at  which  they  diffuse  from  the  catalyst  particles  to 
the  interstitial  fluid,  were  considered 


G dY . . 

( — ) __li  = k3(Y.  - Y.  .)  A 

p.  dz  c js  jl  p 


j = h , a , Hj , Nj 


Where  G is  the  mass  flow  rate  per  unit  cross  section  of  the 
chamber,  p is  the  density,  Y the  mass  fraction  of  the  different 


species,  z the  axial  coordinate  along  the  chamber,  a mass 


transfer  coeficient  that  depends  on  the  interstitial  fluid 


properties,  and  A^  the  external  surface  of  catalyst  particle 


per  unit  volume  of  bed.  Subscripts  and  superscrips  h and  a 


refer  to  hydrazine  and  ammonia  respectively,  and  i and  s to 


interstitial  fluid  and  particle  surface  conditions. 


A similar  equation  is  given  for  the  energy  conservation 


dT. 

c 6 -r-±-  = h (T  - T.  )A' 
p dz  c s 1 p 


where  c^  is  an  average  specific  heat  at  constant  pressure  of  the 


:gas  mixture,  T the  temperature,  and  h a heat  transfer  coefficient.. 

c 


This  is  a system  of  five-  equations  to  determine  the  temperature 
and  mass  fractions  of  the  four  species  in  the  interstitial  fluid 
along  the  chamber.  In  order  to  complete  the  system  an  equation 
of  state  is  required,  that  we  take  to  be  of  the  form 


P = pTR 


J8  17  32 


where  the  pressure  P is  assumed  to  be  constant,  that  means 
that  we  neglect  the  pressure  variations  compared  to  the  absolute 


value  of  the  pressure  ^ ^ ; and  also  are  required  relations 


between  the  particle  surface  properties  and  the  interstitial 
fluid  properties  that  we  are  going  now  to  present. 


In  ref.  2 an  analytical  solution  was  found  of  the  problem 
of  calculating  the  temperature  and  mass  fractions  distribution 
inside  the  particle  . by  solving  the  heat  and  mass  conservation 
equations  of  the  four  species, subject  to  the  conditions 
expressing  regularity  at  the  particle  center,  and  that  the  heat 
and  species  fluxes  at  the  su?face  are  those  given  by  the  right 


,(*)  The  validity  of  this  assumption  was  checked  in  ref.  2. 
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hand  side  of  equations  (5)  and  06).  From  this  solution  the  re- 
quired relations  between  the  interstitial  fluid  properties  and 
the  particle^surface  properties  are  found.  The  details  of  this 
resolution  can  be  found  in  ref.  2,  here  we  shall  only  reproduce 
the  final  results.  In  ref.  1 this  particle  problem  was  solved 
numerically,  simultaneously  with  equations  (5)  and  (6)  along  the 
chamber.  In  figure  2 is  presented  an  schematic  of  the  temperature 
and  hydrazine  and  ammonia  mass  fractions  distributions  inside  the 
catalyst  as  they  where  calculated  in  ref.  2.  One  of  the  main 
results,  also  appearing  in  ref.  1,  is  that  the  hydrazine  .mass 
fraction  at  the  particle  surface  is  very  small  compared  to  the 
corresponding  interstitial  value;  the  decomposition  rate  is 
controlled  by  external  diffusion  to  the  particle  because  the 
velocity  of  reaction  of  hydrazine  is  very  large.  Consequently, 
a detailed  knowledge  of  the  kinetic  mechanism  of  hydrazine 
decomposition  is  not  required  for  this  study.  Whereas  hydrazine 
decomposes  in  a very  thin  region  close  to  the  particle  surface, 
ammonia  decomposes  in  a much  thicker  region  that  may  or  may  not 
be  small  compared  to  the  particle  radius.  The  following  equations 
exprexing  an  species  and  energy  balance  at  the  particle  surface 
give  the  required  relations  between  the  particle  surface  properties 
®.n4  the  interstitial  fluid  properties 


p.KJ(Y..  - Y.) = - v . 
i c 3i  ]S  aj 


p D Y 
s as  as 


A-v.  .p  .K  Y.  . 
h]  i c hi 


h (T.  - T ) = H P-s  as-as.  A - H.Khp.Y  . 
ci  s a a hcihi 


’hs  ■ 0 


(10) 


where 


j = a,  N2,  H2 


i 

I 

\ 


- w-— wMrtfffl&& 


and  v^j  are  the  stoichiometric  coefficients  of  the  decomposi- 
tion reactions,  v equals  -1.  The  left  hand  sides  of  (8)  and 

aa 

(9)  represent  the  external  fluxes  to  the  particle  (see  equations 
(5)  and  (6));  the  first  and  second  terms  of  the  right  hand  sides 
represent  the  ammonia  and  hydrazine  decomposition  contributions 
to  the  internal  fluxes,  respectively,  a is  the  particle  radius, 

A is  a non-dimensional  number  proportional  to  the  slope  of  the 
ammonia  mass  fraction  at  the  inner  edge  of  the  hydrazine  decomposi- 
tion region,  and  D is  the  ammonia  diffusion  coefficient  in  the 

3 

particle.  For  an  ammonia  decomposition  rate  of  the  form 


k = b exp(-  T /T) 
a a 


(11) 


where  b is  constant'  , A is  found  to  be  a function  of  the  three 
non-dimensional  parameters 


Y = T /T 
a s 


(12) 


2 2 
# = k a /D  • 

as  as 


(13) 


8 = p D H Y /(K  T ) 
s as  a as  s s 


(14) 


where  K is  the  heat  conduction  coefficient  inside  the  particle. 

The  parameter  8 expresses  the  ratio  of  the  mass  and  heat 
diffusion  velocities  inside  the  particle  and  is  very  small  because, 
whereas  the  diffusion  is  controlled  by  the  gas  inside  the  particle, 
the  heat  transfer  is  controlled  by  the  solid  part,  that  usually 


'according  to  ref.  1,  b depends  on  the  hydrogen  mass  fraction 
like  Yjj^*  . However  we  may  assume  that  hydrogen  diffuses  so  fast 
in  the  particle  that  its  mass  fraction  is  constant  equal  to  its 
surface  value. 


LtUL 


10.- 


is  made  of  aluminum  or  other  good  heat  conducting  metals 
Y3  = 0,  A takis  the  form 


For 


A = 4>  cth  4>  - 1 


(15) 


For  large  values  of  $ and  moderately  small  values  of  y3 


2„2. 


A - *(1  + y3/6  )-( 1 + y3/9)+  0Cl/*r)  + 0 ( Y "3  )+  0(3) 


(16) 


As  we  shall  see  later,  the  contribution  of  ammonia  decomposition 

is  only*  important  for  large  values  of  T , when  § is  large,  and 

s 

we  shall  only  retain  the  first  term  of  the  right  hand  side  of  (16) 


We  now  need  boundary  conditions  for  equations  (5)  and  (6) 
that  give  the  temperature  and  mass  fractions  at  the  beginning 
of  the  vapor  region.  The  initial  temperature  is  the  boiling 
temperature  of  hydrazine  at  the  operating  pressure 


T = Tb(P) 


(17) 


The  initial  value  of  the  hydrazine  mass  fraction, 


is 


determined  by  the  fact  that  enough  hydrazine  should  decompose 
to  liberate  the  heat  necessary  for  the  hydrazine  to  vaporize  and 
reach  its  boiling  temperature.  In  ref.  2 this  initial  value 
is  given  in  the  form 


H1 

Tf  = V+  F 


hb 


(18) 


where  T^  can  be  interpreted,  as  we  shall  see  later,  as  the 
maximum  temperature  that  would  be  reached  by  the  mixture  if 
there  would  not  be  ammonia  decomposition.  Ref.  2 gives  for  a 


= 1 


Tf  = 1450°K 


(19) 


11.- 


Ref.  1 gives,  also  for  o = 1,  the  expresion 


Tf  = 1420(1  + 0 . 035( P/70 ) ) °K 


(20) 


, with  p in  atmospheres,  that  gives  a slight  variation  of 
with  pressure  from  1421°K  for  p = 1 atm  to  1463°K  for  p = 70  atm. 
Refs.  3 arid  4 give  a value  of  1650°K  for  a = 2/3.  Using  the  same 
procedure  of  ref.  2,  and  tdang  for  the  expression  given  in  (3) 
the  following  correlation  seems  to  be  appropiate  fo.r  T- 


T.f  = 1450  + 600  (l-ct) 


(21) 


that  reproduces  the  results  of  refs.  1 to  4. 

The  initial  values  of  the  mass  fractions  of  the  other 
species  are  obtained  as  functions  of  the  initial  hydrazine  mass 
fraction  by  considering  that  no-ammonia  has  been  decomposed  in  the 
initial  vaporization  region.  From  the  stoichiometry  of  equation 
(1)  it  is  obtained: 


Yab  - (2  ~ - Yhb> 


(22) 


YN2b  16  ° (1  Yhb* 


(23) 


fH2b  " 52  (3°  - 2)(1  " Yhb> 


(24) 


Before  proceeding  further, it  should  be  noticed  that  three 
first  integrals  of  the  system  of  five  differential  equations  (5) 
and  (6)  can  be  obtained  by  linear  combination  of  these  equations. 
On  using  the  previous  initial  conditions  (17)  (18)  and  (22)  to 
"(24)  it  is  obtained  for  these  integrals 


■ 


ft. 
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2.  SOLUTION  WITH  NO-AMMONIA  DECOMPOSITION 


•As  It  was  stated  in  the  introduction , the  main  simplification 
of  this  problem  consists  in  considering  that  the  ammonia  decompo- 
sition is  only  important  at  high  particle  surface  temperatures, 
close  the  maximum  one  reached  with  no  ammonia  decomposition. 


Consequently  the  solution  obtained  in  this  part  will  reproduce 


many  of  the  as  poets  of  t/he  general  solution  obtained  in  the  next 


part,  and  w.. . x also  give  the  optimum  value  of  the  specific 


impulse . 


To  neglect  ammonria  decomposition  corresponds  to  make  the 
fir.st  terras  on  the  right  nand  side  of  equations  (8)  and  (9) 
equal  to  zero.  With  this  approximation,. and  considering  equation 
(i!0) , equation  (5)(for  hydrazine) and  equation  (6)  can  be  written 
in  the  form  f 


dY 


hi  _ 


dz 


- K-1  A 
c p 


lhi 


(28) 


dT : 


c G ~ 
p dz 


= - H.K'p.Y.  . 
h c i hi 


(29) 


By  dividing  these  two  equations  it  is  obtained 


dT , 


c = - H. 

p dY*,  • n 


(30) 


hi 


and  on  integration  with  the  initial  condition  (18) 
H. 


Ti  = Tf  - 


hi 


'When  all  hydrazine  is  decomposed  T^  reaches  T^. 
-goes  to  zero  exponentially  like 


(31) 

From  (28)  Y 


hi 


~~—&£  : -tr-  i *'  1 j - |~|-  1 I “*  ^ **  1 •■** 


.31 


“S 


3 

•£ 


I 
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3 

% 


£ 

4 

4 


Yv . = Y..e 
hi  hb 


(32) 


where 


« = /:  kc  pi  dz 


(33) 


Kp  . is  an  slowly  varying  function  of  T.  that  could  i..  principle 
be  considered  constant  with  a certain' average  value  of  between 
and  Tf  . 

By  taking  (31)  into  equations  (25)  to  (27), we  obtain  the 

concentration  of  the  species  as  functions  of  Y,  . 

ni 


Yai  * M <2-.)(t.TM) 


(34) 


= it  °(1-Yhi> 


(35) 


yh2.  = H r (1-Yhi> 


(36) 


For  Y^  equal  to  zero  we  obtain  the  conditions  when  all  j-ydrazine 
is  decomposed 


Tf  * Yaf  S 32  * ^N2f?  16  ° ’*  ^H2f 


3a-2 

32 


(37) 


The  corresponding  specific  impulse  when  this  gas  is  expanded  to 
the  vacuum  is 


2 _ 2hf  _ 2 r . . v 
lf  ’ ■ 2 - ~2  1 h3f  Y3f 
g . g 


(38) 


Using  appropiate  interpolation  formulas  for  the  enthalpies  of  the 
different  qpscies  taken  from  ref.  ,6 


= (56.49+2.44a)Tf  - (20117  - 8125a) 


(39) 


m 


and  taking  for/Tf  the  expression  (21) 


I./.C27.0  see.)  = (l  + 0.326(l-u)  - 0.02(l-<r)*)  * 


* 1 + 0.163(1  - a) 


(40) 


The  values  of  the  particle  surface  properties  are  obtained 
from  equations  (8)  and  (9),  with  the  ammonia  production  terms 
equal  to  zero,  and  the  values  of  the  interstitial  properties 
calculated  before  in  (31)  and  (34)  to  (36) 


Y = 17/32  (2 

ao 


Y„  = 7/16 
N2s 


- 3a-2 


EK 

1.t  ( f- - i)Yhi 

- c - 

“ F + ( K^“  " 15  Ym] 
[*.♦  <K%  - i)Yhl] 


(41) 


(42) 


(43) 


H.  p . K c 

T = T 1 + — - — - ( ■ c -E 

As  Af  . c T-  ' h 

p f c 


- Y>Yhi 


(44  ) 


For  all  K equal  and 
c 


(45) 


, as  it  was  assumed  in  ref.  2,  the  surface  properties  do  not 
change.  If  Le  > 1 Tg  starts  with  a value  higher  than  T^. 
However,  in  general  Le  is  smaller  than  one,  around  0.6. 


T SkV  SRv.  , 


- 16,- 


1 the  hydrazine  is  decomposed  the  particle  surface 
d the  interstitial  properties  become  equal. 
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3..  SOLUTION  WITH  AMMONIA  DECOMPOSITION. 

The  ammonia  decomposition  effect  is  represented  by  the 

first  term  oh  the  right  hand  side  of  equations  (8)  and  (9).  If 

we  make  equation  (9)  non-dimensional  by  dividing  it  by  h Tf,  the 

c r 

relative  importance  of  the  ammonia  decomposition  is  measured  by 
the  term 


i p D Y 
a rs  as  as 

h a T* 
c f 


The  first  factor  in  backets  is  snail,  of  the  order  of  10~  , so 
if  this  effect  is  going  to  be  significant,  A has  to  be  large 
and  the  parameter  $ must  be  large  too.  According  to  this-,  we 
take  for  A the  expression  given  in  equation  (16)  keeping  only 
the  first  term  of  the  right  hand  side.  On  using  equations  (11) 
and  (13),  we  can  rewrite  the  ammonia  decomposition  term  as. 


J exp  ( 


T -Tc 
s f 


2T  T 
f s 


where,  for  typical  values  of  catalytic  reactors, 

1/2  1/2  T 

J.=  Ha  Os  »as  Yas  bs  **>(  ' 


(46) 


T-  T -Tj. 

• g s f 

is  of  order  one.  On  the  other  hand  the  factor  exp(^= ^ — ) 

f s 

experiences  changes  of  order  one  for  relatively  small  changes  of 

T such  that 
8 


Tf-Ts  ■'  If 

T T 

8 a 


<<  1 


since  T is  much  larger  than  T,.  According  to  Kesten1  T has  a 

3 a 
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value  of  27778°K.  On  the:  other  hand, if  T-  is  smaller  than  T,. 

'■•7 - - - , - s ' - r 

.by  an  order  of  magnitude  .of  one, the  ammonia  production  term 
becomes  very  small.  In  the  previous  part  it  was  seen  that  with 
n'd-ammonia  decompositiph  the  particle  surface  temperature  was 
given1  by 


H. 


T-  = T (1  - — £ 
s - 1 fV1  cT 


(1  - Le  ) Y‘  . ) 


Pf 


hi 


If  Le  < 1,  Ts  will  be  smaller  than  Tf  in  a quantity  of  order 


one. ( if.  Y^.£.  is  also  of  order  one).,  and.  it  will  be  consistent  with 


the  analysis  to  neglect  the  ammonia  decomposition  effect . However,, 
as  becomes  smaller,  T-  gets  closer  to  Tf , and  the  ammonia 
decomposition  effects  should  be  taken  into  account;.  Consequently 
we  can  assume  that  the  ammonia  production  term  is  only  important 

for  values  of  T close  to  Tl,  and  if  J does  not  exhibit  such  an 

• • s - x 


strong  dependence  on  T as  the  exponential  term,  we  can  make 

” * ...  "5, 


the  approximation  of  taking  for  J,  given  in  equation  (46),  the: 
value  J^.,  corresponding  to  take  for  the  variables  determining  J 
the-  values  given  in  equation  (37).  Equation  (9)  is  then  written 
in  the'  form 


T.-T 


T'-T, 


H. 


= Jf  expv-2T^ 


-)- 


L„Y,  . 


c'Tl  e hi 
P f 


(47) 


This  equation  has  to  be  solved  together  with  those  determining  ! 

the  evolution  of  the  hydrazine  mass  fraction,  and  the  interstitial 
temperature 


dT-i 


Le  tr  = Ts  - Ti 


(48) 

(49) 
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where  the  nony'dimensional  distance  £ is  defined  in  equation  (33) . 


By  eliminating  between  equations  (47)  and  (49)  and  taking  Y^ 
as  independei/t  variable,  a single  differential  equation  is  obtained 


for  T 


s 


lhi 


^ 1 


+ J 


f 2T , 


exp( 


Ta  VTf 


2T.  Tf 


:Ki 


exp 


T 

4 

Jt, 


T -T  - 
s . f 


Ha 


fj  V f 


c;T,  (l-L^Yhi 


(50) 


to  be  solved  with  the  initial  condition 


Y,  . = Y...  , 

hi  hib  * 


Ha 


Ts  = Tf(l-(1-Le)  Yhb) 


(51) 


'Pi 


This  initial  condition  is  obtained  from  equation  (47)  by 
substituting  in  it  the  initial  values  of  Y^  and  (equs . ( 17 )and 
(18))  and  neglecting  the  effect  of  ammonia  decomposition,  that 
according  to  this  theory  is  exponentially  small  initially. 


In  order  to  reduce  the  number  of  parameters  we  rewrite 
equation  (50)  and  boundary  condition  (51)  in  the  non-dimensional 
variables 


♦ = 


VTf 


2T  T 

4if  if 


(52) 


y = y^ti-i'e  ) 


Ha 


h 

e } c_T 


/(2Tf/Ta) 


(53) 
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-a 


1 


i 


% 


jis.4 Jxt-.zLJt k>.v:  -■ 
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obtaining  thej  equation 


(1,/exp*)|i,i-  al.* 


and:  boundary  condition 
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(54) 


y "=•  1/5 


♦ = -1/d 


(55) 


where 


A = Jf/(2Tf/Ta) 


(56) 


2Tf  H. 

« = ^ /(1-Le)  Y, 

Ta  C?Tf  hb 


(57) 


6 is  an  small  number,  and  at  the  entrance  of  the  reactor  ^ is  a 
large  negative  number  that  makes  the  term  exp  $ in  (54) 
exponentially  small.  The  solution  of  (54)  and  (55)  is  initially 
of  the  form 

♦ = -y 

that  coincides  with  the  solution  given  by  equation  (44)  for  no- 
ammonia decomposition.  Consequently,  boundary  condition  (55) 
will  be  substituted  by 


♦ +y  = Q , y » 


(55  a) 


eliminating  the  parameter  6 from  this  problem.  In  figure  3 is 
presented  the  solution  of  (54)  and  .(55  a)  for  different  values 
of  A and  Le  = 0.6,  giving  the  nondimensional  particle  temperature 


In  figure  3 it  can  be  seen  that  the  particle  surface  temperature 
is  initially  smaller  than  the  interstitial  temperature  and  both 
increase,  reaching  a maximum  first  the  particle  temperature,  and 
at  a later  stage  the  interstitial  temperature.  When  the  maximum 
of  the  interstitial  temperature  is  reached,  both  the  interstitial 
and  particle  surface  temperature  are  equal.  In  ref.  2 it  was 
assumed  that  Le  - 1 , and  the  particle  temperature  always 
eshibited  a decreasing  behaviour,  the  same  is  expected  for  Le  >1 . 
However,  it  should  be  remarked  that  for  Le  >'  1 the  particle 
surface  temperature  starts  with  values  close  to  and  the 
ammonia  decomposition  is  important  from  the  beginning. 

In  figure  3 is  also  presented  the  asymptotic  solution  of 
(54)  for  large  values  of  5.  The  corresponding  asymptotic  form 
of  (54)  is 
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i 

that  corresponds  to  neglect  the  effect  of  hydrazine  decomposition 
compared  to  that  of  ammonia  decomposition  that  is  also  small 


y <<  exp  4 1 


(62) 


The  solution  of  (60)  is 


+ 


(63) 


where  c is  a constant  of  integration  that  has  been  calculated  from 
the  numerical  solutions  and  varies  from  -1.7  for  X = 1 to  2.8 
for  A = 100.  The  solution  presented  in  figure  3 corresponds  to 
’e  =0  that  is  consistent  with  -considering  large.  In  figure  3 
it  can  be  seen  that  if  the  chamber  length  is  moderately  larger 
than  the  distnce  corresponding  to  the  attainment  of  the  maximum 
interstitial  temperature,  the  asymptotic  solution  gives  a good 
approximation  for  both  the  interstitial  and  particle  temperatures 
at  the  chamber  exitv  From  equation  (58)  and  condition  (62)  it 
can  be  seen  that  both  temperatures  are  practically  equal  at  the 
chamber  exit. 


An  important  parameter  that  is  usually  found  in  the  literature 
is  the  fractional  ammonia  decomposition,  (1-F) 


F = Y ./Y  * 
ax  af 


If  rai/(2-“> 


(64) 


Expresing  Y . as  a function  of  T.  «T,  with  equation  (27),  assuming 

that  Y.  . * 0,  and  then  expressing  T as  function  of  £,  through 
nx  s 

equations  (52)  and  (6.3)  we  arrive  ct  the  asymptotic  form  of  F 


F 


“a 


In 


(65) 
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If  we  observe  that  the  quantity 


• C Xc 

e = JLl 

e H 
. a 


2Tf  1 

T TrT 
a ar 


is  small , expression  (65)  can  be  approximately  rewritten  as 

i-e 


F « 


it.  n 
le  L 


(66) 


For  a = 1 and  Shell  405,  e has  a value  of  0.26.  A correlation 
for  F has  been  obtained  in  ref.  1 that  has  many  of  the  characte- 
ristics of  equation  (66),  in  particular,  if  we  take  for  e = 0.28 
the  dependence  of  F on  chamber  length  and  mass  flow  rate  turns 
out  to  be  the  same. 


To  complete  this  work  we  shall  now  give  an  expression, 
relating  the  specific  impulse  that  would  be  obtained  expanding 
to  the  vacuum  as  a function  of  F. 


(Ig)2  = 2ZYjhj  j = a,  N2,  H2  (67) 

where  it  has  been  assumed  that  the  hydrazine  concentration  is  zero. 
With  this  assumption  the  mass  fractions  of  the  hydrogen  and 
nitrogen  and  the  interstitial  temperature  can  be  expressed  as 
functions  of  Y . by  equations  (25)  to  (27),  and  on  using  appropriat 
interpolation  formulas  for  the  enthalpies  of  the  species,  taken 
from  ref.  6,  the  following  expression,  is  obtained  for  the  specific 
impulse  as  a function  of  F (equation  (64))  and  a 

(I/If)2  = 1 - 326 (i"a)  <l“F)[I"0*07(l-«)  + 0.07(2-a)(l-F2 


where  I ^ is  the  maximum  specific  impulse  with  no  ammonia  decom 
position  obtained  in  the  previous  section  " 

If  = 270(1  + 0.163(l-ct))  sec. 

and  F is  given  in  formulas  (65)  and  (66). 
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